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Abstract 

The law of transformation of affine connection for n-dimensional 
manifolds as the system of nonlinear equations on local coordinates of 
manifold is considered. The extension of the Darboux-Lame system of 
equations to the spaces of constant negative curvature is demonstrated. 
Geodesic deviation equation as well as the equations of geodesies are 
presented in the form of the matrix Darboux-Lame system of equations. 

1 Introduction 

It is well-known that the integrable equations have numerous applications in 
geometry. The Korteveg-de Vries equation, sin-Gordon equation, Tzitzeika 
equation, Kadomtzev-Petviashvili equation, Zakharov-Manakov system of 
equations and others are the most famous examples of such type of equa- 
tions. This can be explained by the fact that the mentioned above equations 
have the Lax pair representation, which is equivalent to the condition of zero 
curvature for suitable connections. From this we infer that the law of trans- 
formation of affine connection is a key to understanding the nature of such 
type of equations and their integrability. 

The matrix Zakharov-Manakov system of equations discovered in context 
of formal generalization of Inverse Scattering Transform Method in multi- 
dimensional case [1] after its geometrical interpretation as the Darboux- 
Lame system [2-6] gives an example of such relation between equations 
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and geometry because it corresponds to the simplest case of the law of 
transformation of the flat affine connection. For the spaces of non-zero 
curvature the law of transformation of connection leads to generalization of 
Zakharov-Manakov and Darboux-Lame systems of equations. 

The problem of integration of equations of geodesies has a great sig- 
nificance in geometry. We show here that in a certain case this problem 
is connected with the matrix Darboux-Lame system of equations and its 
generalizations . 

The equation of geodesic deviation is also important in geometry, in 
the theory of Riemann spaces. For surfaces it coincides, in essence, with 
the Gauss equation in geodesic coordinates while in the three- and higher- 
dimensional Riemann spaces it carries essential information about these 
spaces. Solutions of the geodesic deviations equation (Jacoby fields) and 
their properties (e.g., the existence of conjugate points) are related to vari- 
ous important characteristics of Riemann spaces. 

In this paper we also consider applications of the matrix Darboux-Lame 
system to the study of the geodesic deviation equation. 



2 Equation of the law of transformation of affine 
connection 

Here we will present for convenience basic formulas, which will be used in 
what follows. Let M n be the space endowed with the affine connection. 
This means that the components of connections T l - k {x l ) and Y l - k {y l ) in two 
various systems of coordinates (x l ) and (y l ) are connected by the relations 

d2 y k _ T i ix sM_ f k (v sM9]r m 



dx^xi tjy ' dx l lmyy 'dx i dxi 

We can consider these relations as the system of partial differential equa- 
tions for n unknown functions y k of variables x k . In general case this is a 
compatible system of equations and its properties depend on the curvature 
tensor 

~dr n dr t 



+ r?r n — r n r; 



OX L OX n 

which has the law of transformation as the four rank tensor 

, _ dx i dy r dy q dy p - s 
kln dy s dx k dx l dx n rqp ' 
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where Ti are matrices with components Y\ k and R\i n are components of 
curvature tensor of the connection T; . 

The system (|l]) for manifolds of dimension three is the main object of 
our consideration. 



3 Triply-orthogonal systems of surfaces, flat coor- 
dinates and the Darboux-Lame system of equa- 
tions 

The above system of equations in three-dimensional case is the system of 
18 equations for three unknown functions u(x, y, z), v(x, y, z), w(x, y, z), de- 
pending on three variables (x,y,z). It contains 18 coefficients of affine con- 
nection of the manifold and their solutions depend on the specific choice of 
coefficients of connection. With the help of coordinates u(x, y, z), v(x, y, z), 
w(x, y, z) the geometry of the 3-dim space is described. Let us consider the 
simplest of examples. 

It corresponds to the flat space, i.e. when the coefficients of connection 
Y l - k = 0, or the tensor of curvature of the connection T l - k is equal to zero: 

The initial system (|l|) takes the form 

8V -r'.(x')^. (3) 



dx l dxi 13 dx l 
For 3-dim space we have the following system of equations 

d 2 u , du 2 ® u -p3 d u d 2u -pi d u r 2 ® u r 3 ® u 



dxdy L2 dx u dy u dz' dxdz Li dx 16 dy Li dz'' 
d 2 u i du 2 9u 3 du d 2 u x du 2 ^ u r 3 & u 



dydz 26 dx 26 dy lA dz' dx 2 ll dx 11 dy 11 dz' 
d 2 u i du 2 9u 3 du d 2 u , du 2 9u 3 du 

W = 22 dx~ + 22 ^y + 22 9? = 33 dx~ + 33 ^J + 33 3? 

and corresponding equations for the coordinates v(x,y,z) and w(x,y,z). 

The system of equations (^) is very well known in classical Differential 
Geometry. This is the Lame system of equations for the triply orthogonal 
curvilinear coordinate systems in a flat Euclidean 3-dim space. At present 
the half of this system containing only mixed derivatives has appeared in 
an non-evident general matrix form in work [1] in context of the multidi- 
mensional generalization of integrable differential equations but without any 
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applications. An explicit form of this system of equations with its clear ge- 
ometrical sense, simplest solutions and applications was presented by the 
author (see [2-6]). 

Here I briefly review some results of this approach. For the Riemann 
3-dim space in orthogonal metric 

ds 2 = A 2 (x, y, z)dx 2 + B 2 (x, y, z)dy 2 + C 2 (x, y, z)dz 2 

the condition 

Rijkl = K {9ik9jl - 9il9jk) 
leads to six equations, three equations of which do not contain parameter k. 



A zy = (C y /C)A Z + (B z /B)A y , 
B ZX = (C X /C)B Z + (A Z /A)B X , 
C xy = (A y /A)C x + (B x /B)C y . 

These equations named the Darboux system can be considered as the 
scalar reduction of the general matrix Zakharov-Manakov system of equa- 
tions [1] and they can be integrated by the Inverse Scattering Method using 
the following linear problem: 



* zy = (C y /C)* z + (B z /B)* v , 
$** = {C x /C)® z + (A z /A)$ x , 
$ xy = (Ay/A)$ x + {B x /B)$ y . 

The partial case of the Darboux system of equations is connected with 
theory of the normal Riemann space. The notion of normal Riemann space 
was introduced by Eisenhart. 

Definition 1 The n- dimensional Riemmanian space with local coordinates 
u % is normal when the conditions on main curvatures K; L is fulfiled 



du l 

and 



du l 

d In gij 2 dKi 



if I, 



du l K\ — Ki du l 
where A and fi are some functions of coordinates u 
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Remark 1 The values K\ , K 2 , ■ ■ ■ , K n are given the name of the princip- 
ial curvatures relatively to some symmetrical tensor bij of n- dimensional 
Riemmanian space M n with metric 

ds 2 = gijdu l du^ 

if they are the roots of algebraic equation 

\bij - Kgij\ = 0. 

The inherent vectors ^ of main directions of the tensor bij from the 
equations 

(b ij -K h9ij )a = o 

are defined. They are orthogonal 

Qij&l = v + 1, 
and satisfy to the conditions 

b*&i=0, p + q. 

According to [7] the system of equations for the principal curvatures in 
the 3-dim case looks as 

(K 2 - K 3 )K lx + 3(K 3 - K x )K 2x + 3(#i - K 2 )K 3x = 0, 

3(K 2 - K 3 )K ly + (K 3 - K x )K 2y + 3{K X - K 2 )K 3y = 0, 
3{K 2 - K 3 )K lz + 3{K 3 - K x )K 2z + (K x - K 2 )K 3x = 0. 
Using the relations 
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we get 



"' 1 K ly =0, 
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or 

By analogic way the system of equations for and K3 have been writ- 
ten. 

For the construction of partial solutions of this system we can use the 
linear system of equations 

+ W^T) I** ~ *»] = °> ( 4 ) 

**» + 2(^J " $ yl = °" 

If tp and x are two solutions of the system (01] connected by the relation: 

&u~ =x/2 ~ Ul du~' (5) 

where u% = (x,y,z), then the equations (Qj are the conditions of compati- 
bility for 

So, if the function Xh is a solution of the system fy), then the function 
(fh also will be a solution. As result we have the following relations between 
the solutions 

dXh+i /0 dxh 

OUi OUi 

As example, begining from the trivial solution xo = we can obtain the 
solutions in the form of symmetrical functions of variables (x, y, z) [5, 7j. 

Recently the Darboux and the full Lame system of equations have been 
integrated by the IST-method [8,9,10]. 

Remark 2 The full Lame system of equations can be applied to integration 
of the Einstein equations. The simplest illustration of that is the result of 
E.Kasner on the representation of the Schwarzschild solution of the Einstein 
equations 

ds 2 = (1 _ ^ )df 2 _ (1 _ ^yl dr 2 _ r 2 d() 2 _ r 2 gin 2 ^ 2 

r r 



6 



as six squares of differentials in 6- dim flat space 



ds 2 = -dx 2 - dy 2 - dz 2 + dX 2 + dY 2 - dZ 2 , 



where X, Y, Z are defined by equations: 




1 + 



(R 2 + 16m 2 ) 3 



256m 4 



dR, X = 



VR 2 + 16m 2 ' 



Rsint 



Y = 



Rcost 




VR 2 + 16m 2 



Using this result we can obtain the Schwarzschild solutions of the Ein- 
stein equations from the Lame system for the flat six- dimensional space. 

It is apparent that arbitrary J^-dim metric can be presented as embedded 
in a flat space of suitable dimension and use then the theory of the Lame 
equations for integration of the Einstein equations. 

Remark 3 The theory of the Lame system of equations is connected with 
the Euler-Picard system of partial differential equations, having important 
applications to number theory. The equations have the form 



where I and n are natural numbers such that < I < n and the integral 
functions 



of variables t = t\...t r and of the family of cycles at on compact Riemann 
surfaces of planar equations of the type 







Y n = (X- 1)X(X - h)(X - t 2 )...(X - t r ) 



are solutions of the Euler-Picard system [11]. 
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4 On equations for the coordinates of distorted 
spaces 

The simplest generalization of the Lame system of equations is connected 
with the manifolds of constant negative curvature. They have the metric 

2 _ dx 2 + dy 2 + dz 2 
as — ~ 
z z 

and for such a metric the matrices of the ChristoffePs symbols are 
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From the fundamental system of equations (1) for the coordinates we get 

(T{ 2 + HhL) Ux + (T 2 2 + — H, 



/-pi , w z\ i /-n3 i 
( r !3 + — )«« + (ri3 + )Uz, 



III, 



w 

IV 7 



u 



•yz 
Vxy 



W 



(r| 3 + -K + <X% + 
(r} 2 + ^fK + (r? 2 + ^K, 
(r| 3 + -K + (r? 3 + -K, 



^ = (rl 3 + + (lis + 

where the coefficients rf 2 , T 2 2 , and T\ 2 can be reduced to zero using the 
special choice of system of coordinates. 

The first three of above equations have the form of the Darboux system 
of equations and the equations for the coordinate w(x, y, z) can be written 
as follows: 



w x w y 



Wxz + 



w 

w x w z 
w 



Ay 

A 



Wxy + — = -TW x + ~^-Wy ~ 
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Cx 



A z 



W yz + 
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Bz 
B 



Wy + 



Cy 
C 



w z 



U X Uy 


+ V X Vy 
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u x u z 


+ v x v z 
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^y^z 


+ VyV Z 
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After the change of variable w = v R we obtain the following compatible 
system of equations 

B C 
R zy = ~77 R y + ~FT R z ~ 2 ( u y u z + v y v z ), 

-Rxj/ = -j^ R X + ~~ 2(u x 1i 2/ + V X Vy). 

The solutions of this system of equations can be used to study problems 
of the theory of 3-dimensional manifolds, the knots theory and so on. 

It is apparent that in an anologous way one can investigate the Einstein 
spaces and other more general affme-connected spaces. 

Remark 4 It can be shown that the fundamental system of equations 
(and (^) can be integrated (using some modification!) with the help of 
representation of the functions u l in the form 

x y z 
u l (x,y,z) = U\x,y,z) exp[- + — - + — j]. 

Corresponding solutions for connections coefficients and coordinates u 1 
of space have important applications to different problems of Geometry. 



5 On geodesies of the space of affine connection 

The problem of integration of equations of geodesies is very important in 
Differential Geometry. The simplest case of such type of equations is con- 
nected with the second order ODE 

y" + ai(x, y)y' 3 + 3a 2 (x, y)y' 2 + 3a 3 (x, y)y' + a 4 (x, y) = 0. 

It turns out that the theory of the matrix Laplace equation is connected 
with this problem. 

The general equations of geodesies of the affine connected spaces with 
coefficients r,-„- are: 

d?x l , dx k dxi 
ls T+ kj ^s~~ds~ = 
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Let us present them in the following form 

d 2 x « / dx\ 2 . dx dy x / ' dy\ 2 x dx dz 3 

~d^ + ill \Ts) +lil2 TsTs +i22 \d~s) + lil3 TsU 



i dy dz J x dz fc dz- 7 



d 2 y ^9 f dx\ 2 dx dy ^ 9 / dy\ 2 „„ 9 dx dz 3 

tI + T i — +2r 2 2 — -f + r^ 2 -f +2r 2 ,— — + 

ds 2 \ds/ iZ ds ds \ds J 13 ds ds 

2 dy dz J 2 

+2r2 ^d7 + rfe ^^7 = ' 

d 2 z* , / dx\ 2 : dx dy , / dy\ 2 , dx dz 3 

■3-5- + r n ( 3- ) + 2r i2-^-# + r22 -f +2r*i 1 -3-- i -+ 

(is 2 \ds/ 1Z ds ds 11 \ds J 13 ds ds 

+2i2 ^dr + ife ^^" U ' 
writing the coordinates x* in the form x* = (x,y,z l ). 

Then we will consider the coordinates 2* as the functions of variables 
x, y. So the following relations are fulfilled: 

dz 1 j j dy 

ds x ds y ds 



dV , fdx\ 2 , dxdy , /dy\ 2 , d 2 x , d 2 y 



and 



ds 2 xx \dsJ xy dsds yy \dsJ x ds 2 y ds 2 ' 

Putting these relations in above formulas one obtains the system of equa- 
tions for the functions z l (x,y) (after making equal to zero the expressions 
at the derivatives (f ) 2 , (§) 2 and (f g). 



z Lr — [rii+2ry^+r^z^^]z^,-i-[r 2 1 -i-2r 2 J -^-i-r| J -2:^^]2:* 2r^^ r l n r^z^zi, 



— [^12 + ^ij^ + r^-^ + ^kj( z x z y + z y z x)\ z x Jr 



+ [^12 + ^1^^ + r2j z x + ~j^\j( z x z y + 

— r* — r* ^ — r* — -r* (y k ^ + r k ^\ 

L lj z y L 2j^x 1 12 2 k 3^- x y y x ' : 
z yy = [r^+2r2j' z :y+rj^zkz;j]z£-|-^^ 

So one obtains the following statement 
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Proposition 1 There is one-to-one correspondence between the second or- 
der ODE 

y" + oi(x, y)y' 3 + 3a 2 (x, y)y' 2 + 3a 3 (x, y)y' + a 4 (x, y) = 0. 

and two-dimensional surfaces of the affine connected space A n . 

The following relations between the coefficients ai(x,y) of the equation, 
the coordinates z k {x,y) of the space A n and the coefficients of the connec- 
tions Tfj(x,y) are true 

y" - [T\ 2 + 2Y\ jZy + T\ jZ k y z y ] y' 3 + [T 2 22 - 2F\ 2 + Y%z k y z y - 21^4 - 2T\ jZy + 

+2T 2 jZ y — T\^{z k z x + z x z k )\ y' + 
+ [2T 2 12 - rh + r%-(44 + ziz k y ) + 2T%zi + 2T%z y - 21^4 - y'+ 

+ 2TfjZ 3 x + Tf.jZ^.zi = 0. 

The above mentioned system of equations for coordinates z l (x,y) in 
general case is nonlinear generalization of the matrix Laplace equations. 

Restricting our consideration to the linear part of this system of equa- 
tions we get 

Z xx = A(x, y)Z x + B(x, y)Z y , 

Z xy = C(x,y)Z x + D(x,y)Z y , (7) 
Z yy = E(x, y)Z x + F(x, y)Z y , 

where Z = z l (x, y) is a vector-function, A, B, C, D, E, F are matrices. 
From the conditions of compatibility the relations follow: 

A y - C x = [C, A] + DC — BE, 

D y -F x = [FD] + EB- CD, (8) 

D x - By = AD + BF - CB - D 2 , 

C y - E x = EA + FC - DE - C 2 . 

Let us consider some particular cases. 
1. N = 3. The linear system is: 

Z X X = ~ 2rf 3 )^a; + T\lZy, 
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z xy = (rj 2 — T?a)Zx + (r?2 — rfg)^, 

■^2/2/ = ^22^x + (F 22 ~~ ^T^ajZy. 

From the condition of compatibility = and Z yyx = we get: 

r?3 = 3(^11 + ^12)) = -(rj 2 + r| 2 ), r 23 = r;[ 3 

and as result we obtain the system of equations 

Z xx = 3T1I ~~ 2r^ 2 )-^a; + ^flZy, 

Z X y = 3(^12 ~~ ^%i)Z x + -(2Tl 2 — ^\l)Z y , 
Zyy = ^\lZx + 3(^22 ~~ 2T\ 2 )Z y . 



The corresponding equations of geodesies are 

y 

From the conditions of compatibility one obtains the equations 



' ^2/ + (rl 2 - 2T\ 2 )y' 2 + (2r? 2 - T^y' + I* = 0. (9) 



Q-ix — &2y — 2a^,a\ — 2a|, a2x — ^3?/ — ai«4 — 02^3) 

0.4y — fl3x = — 2a2d4 + 2a 3 . 

These conditions for coefficients cij correspond to the ODE (9) with the 
projective flat of connection, i.e. the components of its curvature tensor are 
equal to zero [12,13]. 

2. N=4 In linear approximation we obtain 

z xx = (Tn — 2rf 3 )z^ + ^\\Zy — 2T\ A z*, 

7 3 _ / r l _ r 3 \ 3 , (V 2 _ r 3 \ 3 _ V 3 4 _ r 3 4 
z xy — \ L 12 i 23/ z i: i ^V i 12 1 13/ z ;y 1 24 z x 1 14 z y 

Z TO = ^22 z x + (^22 ~~ ^23) z y — 2T\^Zy^ 

z ix = (Tii — 2rf 4 )z^ + T^Zy — 2r^ 3 z^, 

A _ (r l r 4 ^4, (r 2 r 4 n 4 r 4 3 p 4 3 
2 a;i/ — I 1 12 — 1 24j z :r + V 1 12 — 1 14j z ;y ~~ 1 13 z ?/ ~~ 1 23 z x; 



z 4 

to 



^22^ + (^22 ~~ 24) Zy — 2r 23 ^. 
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Or in the matrix form 

/ rh - 2r? 3 -2r? 4 \ / r 2 u o , 



-2T\ Z T\ 1 -2T\ i r^\ VI, 



7 _ I r 12 r 23 r 24 1 7 I [ r 12 r 13 r 14 , 

^2/ — p4 pi "p4 I 3* ' I t4 -n2 -n4 ^'y 

\ _i 23 1 12 ~~ 1 24 / \ — 1 13 1 12 ~~ 1 14 ' 



7 [ r 22 o \ / r| 2 - 2r| 3 -2r| 4 \ 

o r^ 2y )^ + l v -21* r| 2 -2r| 4y )^- 

Let us consider the surfaces in 4-dim space corresponding to the equa- 
tions y" = f(x,y). According to (9) one get 

l/' + Tii = 0, 

and the following relations are true: 

Tl 2 = 0, T 2 22 - 2T\ 2 = 0, 21^-1^ = 0, r 2 u ^o. 

From these relations we get 

F = 2C, A = 2D, E = 0, 

and the system takes the form 

C x - 2D y + 2CD - DC = 0, D x - B y + CB - 2BC - D 2 = 0, 

D y - 2C X + 2DC -CD = 0, C y = C 2 . 
So from the first and third equations we have 

C x - D y + CD - DC = 0, 

and the matrices C and D can be presented in the form 

C = Q y Q~ 1 , D = e x o~ 1 , 

where @(x,y) is matrix satisfying the compatible system of equations 

e X y = e x Q- 1 e y + QyQ- 1 e x 

Qyy = 2QyQ- 1 Qy. 
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The foregoing method of investigation of equations of geodesies can be 
generalized for the spaces of higher dimensions. So in three-dimensional case 
we obtain the following equation of geodesies 

y" - Tl 2 y' 3 + (T 2 22 - 2T\ 2 )y> 2 + (2T 2 12 - T^y' + T 2 n + 

+2IV + IV 2 + 2(T 2 23 - r\ 3 )y'z' - 2T 1 23 y' 2 Z ' - T l 33 y' z' 2 = 0, 

z" - Tl 3 z' 3 + (r| 3 - 2T\ 3 )z' 2 + (2r? 3 - r^y + r?i+ 

+2T 3 12 y' + T 3 22 y' 2 + 2(r| 3 - F\ 2 )y'z' - 2V\ 3 y' z' 2 - T l 22 y' 2 z' = 0. 

Prom these relations the matrix Darboux-Lame system has appeared by 
natural way. 

In fact, the coordinates x % can be written in the form x,y, z,q l (x,y, z). 
Then the following relations are fulfilled: 

< ^=q i = q i x x + q y y + q\z 

and 

d 2 q i 

Putting these relations in the general equations of geodesies: 

(fx 1 , dx k dx j _ 

we obtain the system of equations for the functions q l (x,y) (after making 
equal to zero the expressions at the derivatives (^f) 2 , (^f) 2 > (^f) 2 ; an d 

V ds" yds ds^ {ds dsr 

It has the form: 

Qlx = [Th + + V\ j q k x q x \q x + 

+{v 2 l + 2T 2 j qi + T 2 j q k x q x ]q y + 

+ [r?i + 2T%qi + \yU- 'l' - 2T\ 3 qi - Y'^U ~ Hi, 
Qly = \Fi2 + Iygj + ^4 + -Vl^qlql + q k q x )]q l x + 
+ [r? 2 + T 2 ljq i + + + q*q> x )]qi+ 
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+ [r? 2 + T\d + r%qi + l -T%{q k x qi + g*g£)]g»- 
4 = t r 23 + ^ + + ^r^-(g*^ + q k z qiM + 

q \ z = [r' 3 + 21%^ + 1'/,;'/'' '/' '/.'/ + [ri 3 + 21%.^ + v\ fl k JM y + 
+[ri 3 + 21* <g + - 21^ - vy^i - r 33 . 

So, in linear approximation we get the matrix Lame system of equations 
on coordinates q l (x,y,z). 

6 Equation of geodesic deviation and the matrix 
Darboux-Lame system 

The general equation of geodesic deviation is 

£>V , dx k dx m ■ 

ds2 t- it fcjm ds ds rr u, 

where i?^ m is the tensor of curvature of manifold. 

We transform this equation into an easy-to-use form. According to the 
rule of covariant differentiation we obtain the relation 

Drf drf , ,dx* , 
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and then 



D 2 rf _ Dp' _ dp' j , dx m 
ds 2 ds ds lmP ds 



d 2 rf <^ % jk j dx k dx l { drf dx k , • d 2 x k { drf , • dx k dx r 



+- 



ds 2 ' dx l '' ds ds ' " jk ds ds " lh 
Furthermore, using the equations of geodesies 



ds ds 



d 2 x k h dx p dx q 

+ r = o 



ds 2 



pq ds ds 



we can get the relation: 
d 2 rf 



ds 2 
dx k dx 



i dx m drj 

+ 2Tlm ^Tls- + 



(ix p dx 9 
HJl— — + R] 



d^)k dx k dx 1 
dx 1 ds ds 

dx k dx m 



Putting here the expression for the curvature tensor 



rf = 0. 



pi 

we find the relation 



or 



km 



dT 



dxi dx r ' 



I pi pn pi pro 

raj fcra ram fci' 



nm fcj' 



d 2 rf , 0F i dx m <V 

+ Zi ZraT 



ds 2 



+ 



or 



kl 



ds ds 



dT 



jk 



dx 1 



i "pri -pi pn I 
+ 1 raZ 1 jfc ~~ 1 in 1 fcZ + 



dxi 



I pn pn 



1 nZ 1 fcj 



cix fc dx' 
ds ds 



ry 



whence after rearragement the desired formula for the equation of geodesic 
deviation follows 



d 2 rf 
1? 



j dx m drf dT l kl dx k dx 1 ■ 

+ lm ^kTTs~ + l)xT^tels~ 7] = 



(10) 



This form of equation of geodesic deviation is convenient in applications 
(see [14]). Let us consider some of them. 
In three-dimensional case one obtains 



drf 
ds 



ff = r^x + rfy + rf z z 
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and 

Putting these relations in (10) and using the equations of geodesies we get 
the matrix Lame system of equations 

rf xx + ir\pi-T\ 1 rf x 
rjyy + 2T 2l rf y - Y\ 2 ri 
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